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Abstract 

We present the first results for the static quark potential and the light hadron 
spectrum using dynamical fermions at /3 = 5.2 using an 0(a) improved Wilson 
fermion action together with the standard Wilson plaquette action for the gauge 
part. Sea quark masses were chosen such that the pseudoscalar-vector mass ratio, 
nips/my, varies from 0.86 to 0.67. Finite-size effects are studied by using three 
different volumes, 8® • 24,12® • 24 and 16® • 24. Comparing our results to previous ones 
obtained using the quenched approximation, we hnd evidence for sea quark effects 
in quantities like the static quark potential and the vector-pseudoscalar hyperhne 
splitting. 



1 Introduction 


Recent years have seen a lot of progress in nnderstanding the spectrnm and the decays of 
hadrons using numerical simulations of lattice QCD (for recent reviews see While 

much effort has been invested in controlling systematic errors such as finite-size effects 
and lattice artefacts, the inclusion of quark loops in the stochastic evaluation of the Feyn¬ 
man path integral still presents a major challenge. Therefore, most simulations rely on 
the quenched approximation for which the systematic errors incurred by neglecting dy¬ 
namical quark effects cannot be assessed. However, recent progress in the development of 
efficient algorithms @] and increased computer power, has greatly increased the prospects 
for simulations with dynamical quarks. Several groups have published results for hadronic 
quantities in the light quark sector and the static quark potential from dynamical simula¬ 
tions i-0. At the same time it has been demonstrated that leading lattice artefacts of 
order a in physical observables can be eliminated through the non-perturbative implemen¬ 
tation of the Symanzik improvement programme |^,|^. The non-perturbatively 0{a) 
improved Wilson action has been determined in the quenched approximation [^,|T^, and 
first results have also been reported for nj = 2 flavours of dynamical quarks |^. This 
enables one to study the effects of the inclusion of dynamical quarks whilst having better 
control over discretisation errors. 


In this work we report on calculations of the light hadron spectrum using two flavours 
of 0(a)-improved dynamical Wilson quarks at /? = 5.2. For the improvement (clover) 
coefficient we have used Csw = 1.76, a preliminary estimate kindly supplied by the Alpha 
Collaboration prior to the final result Cgw = 2.017 at /? = 5.2, presented in |^. The 
main aim of this study is to understand qualitative features of dynamical simulations 
with improved fermions, such as the sea quark mass dependence of observables, finite- 
volume effects and estimates of autocorrelation times. The complete removal of 0{a) 
lattice artefacts is therefore not our highest priority in this work. Results for the hadron 
spectrum obtained with the “correct” values of Cgw will be published elsewhere [IK . 


The plan for the remainder of this paper is as follows. In section we describe the details 
of our simulation, including the implementation of the algorithm and the analysis of 
autocorrelations. Section ^ contains our results for the static quark potential. The results 
for the light hadron spectrum and discussions of finite-volume and sea quark effects are 
presented in section Section ^ contains our conclusions. Finally, in Appendix ^ we list 
our results for hadron masses on all lattice sizes and parameter values. 


2 The simulation 


In this section we fix our notation, describe the details of the implementation of the 
Generalised Hybrid Monte Carlo (GHMC) algorithm |l^ on the Cray T3E, and give an 
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overview of the simulation parameters used in our calculation. We end this section with 
a discussion of autocorrelations. 


2.1 Lattice action 


The lattice action can be split into a pure gauge part Sg and a fermionic part Sp 


where 


= Sg[U] + V’], 




is the Wilson plaquette action, and Sp is dehned by 


Sp[U,^p,^jJ] = S^[U,^jJ,^jJ] 'tp{x)(T^,yF^,y{x)'4){x 


( 1 ) 

( 2 ) 

( 3 ) 




Here, is the standard Wilson action, and Cs„ denotes the improvement coefficient 


multiplying the Sheikholeslami-Wohlert term The bare parameters of the theory are 
the gauge coupling [3 = &/and the hopping parameter n. Here we work with a doublet 
of degenerate dynamical Wilson quarks, and hence gauge conhgurations are characterised 
by the set of parameters (/3, k; Csw)- For a description of the GHMC algorithm it is useful 
to rewrite the fermionic part of the action in terms of a complex, bosonic pseudofermion 
field 0. In matrix notation we have 


Sp = (j)i 0 — 2 ^ Indet 


even sites 

The odd-even preconditioned fermion matrix M is given by 

M — 4 — n A~^n 

ivixy -^xx -^xz-^zz zyi 


( 4 ) 


( 5 ) 


where D is the Wilson-Dirac matrix and A denotes the matrix for the Sheikholeslami- 
Wohlert term. 


2.2 Implementation of the GHMC algorithm 

The main limitation of the performance of the Cray T3E is the memory bandwidth. The 
increased complexity of the memory system and the fact that the processors support 
multiple instruction issue leads to a loss of performance even for highly optimised Fortran 
code. We have therefore chosen to write key routines in assembler, whilst using Fortran 90 
for less CPU-intensive parts. Less than 25% of the required run time is spent executing 
Fortran code. 
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Using 32-bit instead of 64-bit precision to represent the fields increases the speed by a 
factor 1.7. However, this gain has to be weighed against the degradation of the acceptance 
rate, reversibility, and the accuracy in the evaluation of the (global) energy difference 
required in the Metropolis accept/reject step. The last issue has been addressed by 
evaluating energy differences site by site and performing the subsequent summation in 
higher precision. 


The inversion of the fermion matrix was performed using the Stabilised Biconjugate Gra¬ 
dient (BiCGStab) algorithm with odd-even preconditioning |^T|. The gain compared to 


using the ordinary Gonjugate Gradient algorithm was 40%. Further algorithmic improve¬ 
ments applied in this simulation are described in [p^]. 


The version of the Gray T3E used in this work consisted of 96 processors, each capable 
of 900 MFlops peak speed. Using all our algorithmic improvements and exploiting the 
architectural features of the T3E we typically achieve sustained speeds of 25-30 GFlops 
on such a configuration. 


We now describe the integration schemes used in the molecular dynamics part of the 
GHMG algorithm. As usual one introduces a set of conjugate momenta P for the gauge 
fields U. The HMG Hamiltonian is then defined as 


n^T{p] + Sclu] + SFlu,4>\4’], 


( 6 ) 


where T[P] is the kinetic energy, and we have written Fi? in terms of pseudofermion helds 
0 and 0h T[P] is related to the evolution operator T in molecular dynamics time r, so 
that for any given setP andP of gauge helds and conjugate momenta, respectively 


/Tp\ /t/(r + dr)\ 

[Tu)' [p{t)) [P{T + dT))' 


(7) 


where dr denotes a hnite interval in simulation time. The operators Tp and Tjj are dehned 
by 


Tp{dT) : U — 

Tp(dr) : P-^P-^dT^{SG + SF). (8) 

Since T represents a numerical integration of the equations of motion it does not conserve 
Ti but introduces an error ATi. For a single application of T this error is expected to 
grow as a power of the timestep dr p4|,p5f| 

AP oc [dry. (9) 

Verihcation of this relation for a given integration scheme provides a check for the correct 

implementation of the equations of motion. We have compared three integration schemes 
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scheme 

(dry 

q 

Ti 

{dr? 

1.982(4) 

T2 

[drf 

3.053(2) 

Ts 

{drf 

5.056(6) 


Table 1: The leading variation of ATi for the three integration schemes considered and 
measnred valnes of q. 


defined by 


Ti = Tp{dT)Tu{dT) (10) 

T2 = TpQ)Tu{dT)Tp{^) ( 11 ) 

Ts = Tp(|dr)T(7(adr)Tp(^^^dr)Tf/(Mr)Tp(^^^dr)T(7(a(ir)Tp(^dr), (12) 

where a = 1/(2 — 2 ^A) and b = —2^/^/(2 — 2^/^). Note that T 2 is the standard leapfrog 
integration scheme. One expects that Ti, T 2 and T 3 canse AH to vary as {dr)"^, (dr)^ 
and {dr)^, respectively. This can be compared to the valnes of q obtained from the 
slope of InATi as a fnnction of ln((ir). Snch a comparison is shown in Tabled The 
nnmerically determined valnes of q agree well with the expected behaviour of AH for the 
three integration schemes, and thus we conclude that the integration of the equations 
of motion has been implemented correctly. In our production runs we have chosen an 
integration scheme which, like T 2 , is exact up to order (dr)^. To this end we dehne the 
operators 


Tcidr) 

: P- 

^P- idr-As^ 

Tpidr) 

: P- 

-^P-^drAs, 


and consider the generalised leapfrog integration scheme 


(13) 


(]t 

T(dT) = P(—) 


r«(|)Tp(^)To(|) 


TAt) 


(14) 


This particular scheme allows for a more efficient evaluation of the derivative of S'g + 5'p 
with respect to the gauge field. It reduces to the simple leapfrog scheme for n = 1. In 
practice, we have used n = 1 or 2 . 

Finally we note that although the Generalised Hybrid Monte Carlo (GHMC) algorithm 
had been implemented, we have used standard Hybrid Monte Carlo (HMC) in all produc¬ 
tion runs. 
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L^-T 

/3 

^SW 

#conf. 

4-sea 

hj 

^val 

8^ ■ 24 

5.2 

1.76 

78 

0.1370 

0.1370 

0.1380 

0.1390 

0.1395 





100 

0.1380 

0.1370 

0.1380 

0.1390 

0.1395 





100 

0.1390 

0.1370 

0.1380 

0.1390 

0.1395 





60 

0.1395 

0.1370 

0.1380 

0.1390 

0.1395 


12 ^ ■ 24 

5.2 

1.76 

151 

0.1370 

0.1370 

0.1380 

0.1390 

0.1395 





151 

0.1380 

0.1370 

0.1380 

0.1390 

0.1395 





151 

0.1390 

0.1370 

0.1380 

0.1390 

0.1395 





121 

0.1395 

0.1370 

0.1380 

0.1390 

0.1395 





98 

0.1398 


0.1380 

0.1390 

0.1395 

0.1398 

16^■24 

5.2 

1.76 

90 

0.1390 



0.1390 

0.1395 

0.1398 




100 

0.1395 



0.1390 

0.1395 

0.1398 




69 

0.1398 



0.1390 

0.1395 

0.1398 


Table 2: Summary of simulation parameters and statistics for the computation of hadronic 
observables. 


2.3 Simulation parameters 


Our simulations have been performed aX (3 = 5.2. In order to be able to study consistently 
hnite-volume effects and the dependence on the dynamical quark mass, we have used 
Csw = 1-76 in all calculations described in this paper. The standard lattice size in this 
simulation was 12^ ■ 24 and was chosen such as to guarantee a spatial volume in physical 
units of more than (1 fm)^. Smaller and larger lattices of size 8^ ■ 24 and 16^ ■ 24 were used 
to monitor hnite-size effects. 


In order to distinguish the bare quark mass used in the generation of dynamical con- 
hgurations from that used to compute quark propagators for hadronic observables, we 
introduce the notation to denote the hopping parameter of the doublet of dynamical 
quarks, whilst reserving k™* for the valence quarks. In order to study the dependence of 
observables on the sea quark mass, gauge conhgurations have been generated at several 
values of In Table ^ we list lattice sizes, the values of and k™* and the number 
of conhgurations. 


sea ^vau combined to 


Quark propagators were calculated for every combination of (k 
form hadronic two-point correlation functions. In order to increase the projection onto the 
ground state, the quark propagators used to form hadronic two-point functions have been 
“fuzzed” at the source and/or sink according to the prescription dehned in [^. Statistical 


errors of observables have been estimated using the bootstrap procedure described in 
using 250 bootstrap samples. 


5 















2.4 Autocorrelations 


The determination of autocorrelation times is important in order to achieve small sta¬ 
tistical correlations among the ensemble of configurations and to eliminate the effects of 
insufficient thermalisation. 


The autocovariance of an observable is defined as 


28 




(15) 


where the subscripts on label the the values obtained on successive configurations. In 
practice, the expectation value (...) is replaced by the ensemble average over a finite 
number of configurations. We define the autocorrelation function of by 


pnit) — rQ(f)/rQ(o), 


(16) 


and from the large-t behaviour of pn one obtains the exponential autocorrelation time 

7 -exp 

PQ{t) oc t ^ oo. (17) 

The slowest mode of p^ is thus characterised by which is relevant for the equilibration 
of the system. By contrast, the integrated autocorrelation time depends on the 
observable and is required for the estimation of the statistical error in fl, once the system 
is in equilibrium. It is defined by 

-1 CO -1 oo 

= 9 = 9 + £ ( 18 ) 

^ -oo ^ t=l 

where the latter equality holds since pn^—t) = pn(t). This definition implies that sta¬ 
tistically independent configurations for quantity are separated by In practice 

one has to truncate the infinite sum in eq. ([TB| ) at some finite value tmax- The resulting, 
so-called cumulative autocorrelation time r™™ 

1 ^meix 

+ E Pa(t) (19) 

^ t=l 

is a good approximation to provided that tmax has been chosen large enough so that 
any further increase does not lead to an increase in In other words, a plot of 

versus fmax should ideally exhibit a plateau for large enough fmax- 

In order to obtain reliable estimates for and r™™, autocorrelations should ideally be 
measured using ensembles containing many more configurations than the value of 
This requirement is not easy to satisfy in simulations whose primary aim is to compute 
hadronic properties, i.e. for which the calculation of observables requires a non-negligible 
amount of CPU time. A convenient quantity to determine autocorrelations is the average 
plaquette, which in our simulations has been measured after every HMC update, and 
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Jj. . ^sea ^conf. 

^cum 

8 ^ ■ 24 0.1370 4900 

0.1380 6700 

0.1390 6600 

0.1395 11800 

> 20 35 

44 +33 43+3 

36 53 + 1 

> 57 85 

123.24 0.1370 6000 

0.1380 6000 

0.1390 5600 

0.1395 5100 

26 +5 29 +2 
35 42 +3 
52 43 

51 51 

163 ■ 24 0.1390 3800 

0.1395 4200 

0.1398 3000 

38 37 + 3 

32 +18 27 +24 

32 +^4 32 + 1 


Table 3: Estimates of autocorrelation times for the average plaquette for (3 = 5.2, Cgw = 
1.76 on several volumes. 

not only on the subset of conhgurations used to compute quark propagators. Although 
depends on the quantity the integrated autocorrelation times estimated from the 
plaquette provide a useful guideline for the computation of hadronic observables. 

Examples of our analysis of autocorrelations are shown in Fig.|l|. Here, the statistical 
errors plotted for p{t) and were estimated using a jackknife procedure. In order to 
take into account the effects of autocorrelations in the error estimate of p{t) and r™™ 
themselves, the original data for these quantities were grouped in bins of size h. Jackknife 
averages were then formed for varying bin size, and by increasing h until the jackknife 
errors stabilised, the error bands in the plots were obtained. Fig.|I|(a) shows a plot of 
In p{t) versus t. and its error are extracted from the linear slope at large t using a 
htting routine. Fig.|^(b) shows the cumulative autocorrelation time. The central value 
of and the error are read off in the region where shows no signihcant variation 
within statistical errors. 

Our results for and estimated from the average plaquette are shown in Table^. 
One observes a pronounced dependence of autocorrelation times with the mass of the sea 
quark. In the range of investigated in our study, and increase by roughly a 
factor of two as one goes to smaller sea quark masses. Also, a mild volume dependence of 
and r™™ is observed, so that autocorrelations appear to be slightly weaker on larger 
lattices. In view of the large errors, however, this dependence is not really signihcant. 

We have calculated quark propagators on conhgurations separated by 60 sweeps on both 
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(a) 



(b) 



Figure 1: Autocorrelations for the average plaquette on 12^ ■ 24 and = 0.1370. (a): 
Inp(t) plotted versus t; (b): the cumulative autocorrelation time plotted against tmax (see 
eq. (|l^). The solid lines follow the computed values for In pit) and The dashed 

lines represent the error bands estimated from a jackknife procedure as described in the 
text. 
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8^ ■ 24 and 12^ • 24, and 40 sweeps on 16^ ■ 24, respectively. 


2.5 Scaling of the HMC algorithm and the value of «:crit 

Here we wish to report briefly on a simple method to obtain an estimate of the criti¬ 


cal valne of the hopping parameter, Kent, based on the scaling behavionr of the HMC 


algorithm with the qnark mass. This method is particnlarly nsefnl becanse it can be 
applied independently of an analysis of spectroscopy data, i.e. withont compnting any 
qnark propagators at all. We stress, however, that it serves only to obtain a preliminary 


estimate of Kent, whose actnal valne has to be extracted from the enrrent qnark mass or 


the qnark mass behavionr of the psendoscalar meson. 

Motivated by the idea that the compnter time reqnired for the generation of a dynamical 
gange conflgnration follows a scaling behavionr near the critical qnark mass, we make the 
ansatz 



( 20 ) 


where Nqg is the nnmber of Conjngate Gradient iterations reqnired to invert the fermionic 
part to some given accnracy, and 5 is a critical exponent. 

If the valne of Kent is known, we expect that IniVcc plotted against ln(l/K— 1/Kcrit) 
shonld be linear with slope 5. Conversely, if Kcrit is not known a priori, we can nse several 
trial valnes for Kent, taking the valne which reprodnees the linear behavionr of IniVcc 
as the preliminary estimate of the trne Kent- Snch an analysis is shown in Fig. ^ for 
K®®*^ = 0.136 on 12^ ■ 24. Here a straight line is obtained between Kent = 0.140 and 0.141. 
This procednre can be optimised by performing a linear fit of InA'^cc- For instance, on 
8^ ■ 24 snch a fit yields Kcrit = 0.14004(4). This is to be compared to the valne obtained 
from the quark mass behaviour of the psendoscalar mass described in section^, which 
gives Kcrit = 0.14047 j*)®, which is reasonably close to the value obtained from the scaling 
analysis of the HMC algorithm. 

The procedure outlined in this subsection has its merits because many inversions are 
performed in a typical simulation, and thus a statistically significant value for Nqq is 
easily obtained. Strictly speaking, one should only consider the first inversion of the 
computation of a new trajectory, since this is the only one guaranteed to be performed 
on a physical configuration (i.e. immediately after the global accept/reject step). 

3 The static quark potential 

In this section we describe the computation of the static quark potential using our dy¬ 
namical configurations. The force between static quarks, calculated from the potential 
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Figure 2: In A^cg plotted versus In (1 /k — l/^crit) for several trial values of Kent for = 
0.136 on 12^ ■ 24. 

serves to determine the lattice scale using the hadronic radius vq Furthermore, we 

study finite-size effects and investigate possible evidence for string breaking. 


3.1 General procedure 


The method to extract the potential V{r) from Wilson loops W{r,t) of area |r| x t 
is standard. We have used the algorithm described in to compute “fuzzed” gauge 
links with a link/staple weighting of 2:1 and between 10 and 20 iterations in the fuzzing 
algorithm. Using two different fuzzing levels, we have constructed a 2 x 2 variational 
basis of Wilson loops and subsequently determined the eigenvalues and -vectors of 
the generalised eigenvalue equation 


i,j,k = 1,2. 


(k) 


( 21 ) 


The eigenvector corresponding to at t = 1, to = 0 was then used to 

project onto the approximate ground state [^, ^. This combination of (t,to) turned 
out to be a compromise between good projection properties and the need to avoid the 
introduction of additional statistical noise. The resulting correlator was then fitted to 
both single and double exponentials for timeslices up to f = 8. As a cross-check we also 
performed exponential fits to the full 2x2 matrix correlator. No significant deviations in 
the fit parameters as a result of different fitting procedures have been observed. 
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3.2 Determination of ro/a on different volumes 


The computed values of the potential V{r) can be used to determine the force F{r) 
between a static quark-antiquark pair separated by a distance r = \r\. As discussed 
in , the force can be matched at a characteristic scale tq to phenomenological potential 
models describing quarkonia. More precisely, ro is defined through the relation 


= 1-65, 


( 22 ) 


which corresponds to tq ~ 0.49 fm. Eq. can thus be used together with lattice data for 
the force to extract ro/a, which then yields a value for the lattice scale in physical units. 
This definition has the advantage that one needs to know the force only at intermediate 
distances. An extrapolation of the force to infinite separation, which is conventionally 
performed to extract the string tension, is thus avoided. Hence, the procedure is well- 
suited to the case of full QCD, for which the concept of a string tension as the limiting 
value of the force appears rather dubious, because the string is expected to break at some 
characteristic distance Vb- 


Our determination of ro/a follows closely the procedures described in |^, 
computed the force F^{rj) for orientations d of Wilson loops according to 


We have 


f’d(n) = \d\-'{V{r)-V{r-d)), 
r] = -|ti|(Gi(r)-Gi(r-d)) 

where GL{r) is the lattice Greens function for one-gluon exchange 

. r d^k cos{k ■ r) 

^ ^ ^ (2n-)3 4 sin(%/2) ‘ 


(23) 

(24) 


(26) 


This definition ensures that F{ri) is a tree-level improved quantity [^. In our study we 
have concentrated on “on-axis” orientations of Wilson loops, i.e. where d = (1,0,0). We 
have obtained estimates of tq in lattice units by a local interpolation of F{rj)rj to the 
point defined in eq. ([^). We emphasise that this procedure does not rely on any model 
assumptions about the r-dependence of the force. 


Systematic errors in ro/a were estimated through variations in the interpolation step 
(e.g. by considering a third data point for F{rj)rj besides those which straddle 1.65), 
and also by using alternative fitting procedures in the extraction of the potential (e.g. 
single or double exponential fits, different fitting intervals). We note that the systematic 
error in ro/a, in particular for smaller values of the sea quark mass, is dominated by the 
uncertainty incurred by considering different points in the interpolation step. A summary 
of our results on all lattices and for all values of is shown in Table The configurations 
on which the potential has been determined were separated by 40 HMC trajectories for 
all lattice sizes and quark masses. 
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Table 4: Results for r^/a for different lattice sizes and quark masses. The lattice spacing 
was obtained using tq = 0.49 fm. The hrst error is the statistical, the second an estimate 
of the systematic error as described in the text. 



Figure 3: The static quark potential for the hrst six on-axis separations r/a= 1,... ,6 
plotted against for different lattice sizes. 
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Figure 4: The hadronic scale r^/a plotted against for different lattice sizes. The 

chiral limit is approximately at the left margin of the figure. Solid error bars represent 
the statistical errors, whereas the dotted error bars denote the linearly added systematic 
error. 

The comparison of results obtained on the 8^ ■ 24 and 12^ ■ 24 lattices shows that there 
are pronounced finite-size effects at = 0.1395, whereas for larger quark masses these 
effects are small. The presence, respectively absence, of hnite-size effects is easily recog¬ 
nised in the values of the potential V{r) itself, as shown in Fig. On the other hand, 
there is remarkable agreement in the data obtained on 16^ ■ 24 and 12^ ■ 24, even at the 
lightest quark mass considered. This is illustrated in Fig. where the results of tq/q on 
all lattices are plotted against I/k®®®'. Our findings can be translated into a bound on 
L/ro, above which hnite-size effects in the static quark potential are largely absent at this 
level of precision. From our results we infer that the bound is 

T/ro > 3.2, (26) 

which, for instance, is still satished for < 0.1380 on 8^ ■ 24. This bound, however, 
should not be generalised prima facie to other quantities, in particular the spectrum 
of hadronic states discussed in section^, for which hnite-volume effects could well be 
different. We will return to this point in subsection |4.3| . 

Fig. ^ shows that the data for tq/ a obtained at the three lightest quark masses on L/a = 12 
and 16 show a linear behaviour. We have therefore attempted a linear extrapolation of 
ro/a to the chiral limit using the data at k®®*^ = 0.1390, 0.1395 and 0.1398 only, despite the 
lack of a theoretical motivation as to why such an ansatz for the quark mass dependence of 
ro/a should be valid. Taking into account only statistical errors in rg/a, the extrapolations 


13 




for L/a = 12,16 yield 


ro/a = 4.10 ±0.06, a = 0.122 ± 0.002 fm 
ro/a = 4.15 ±0.06, a = 0.121 ± 0.002 fm. 


(27) 

(28) 


Lja = 12 : 

L/a = 16 : 


Thus, the overall box sizes in the chiral limit amount to 1.46(2) fm for L/a = 12 and 
1.93(3) fm for L/a = 16. A comparison with data for ro/a obtained in the quenched 
approximation shows that our estimates in the chiral limit for nj = 2 massless quarks 
at /3 = 5.2, Csw = 1-76 roughly correspond to values around (3 = 5.85 — 5.9 in quenched 


QCD |3|,|6 


3.3 Is there evidence for dynamical quark effects? 


We now examine our data for the static quark potential for possible evidence for the 
effects of dynamical quarks. In Fig.|^ we show the potential in units of tq, normalised to 
l/(ro), for hve values of used on 12^ • 24. We compare our results to the expression 


[Vir) - l"(ro)]ro = (1.65 - e) - l) - e - l) , (29) 

which follows from the standard linear-plus-Coulomb ansatz for V{r), viz. 


V (r) = Vo + ar — - 
r 


(30) 


and the condition eq. (|^) . Here, a denotes the string tension, and we have set e = 7r/12 


37| , so that the solid line in the hgure has not been obtained through a £t. The data at 
different k®®®' have been offset by H(ro), whose value was obtained by a local interpolation 
of the potential. 


In the presence of dynamical quarks one expects a deviation of the data at large separa¬ 
tions from the linear behaviour described by the curve in eq. (|^) , so that the potential 
in full QCD flattens out due to string breaking. For separations larger than the breaking 
scale rb, one expects that the potential is equal to the mass of two “mesons”, consisting 
of a static quark and a light antiquark, i.e. the energy of a state corresponding to a 
broken string. The masses of such static-light mesons have been calculated on 12^ • 24 for 
^sea _ Q_x 3 gQ g^j^d 0.1395 using the technique described in [^. The error bands of this 
determination are shown as the dotted (k®®®" = 0.1395) and dashed (k®®® = 0.1390) lines 
in Fig.|^. 

Our data for the potential for distances r > 1 fm are neither in disagreement with the curve 
in eq. (p9|), nor with the expected asymptotic value of 2M^*®*'®. In order to check whether 
more statistics could help in revealing the flattening of the potential at large distances, 
we have computed the potential on 194 stored HMC trajectories for k®®® = 0.1395 on our 
larger lattice size of 16^ ■ 24, but no qualitative change compared to the data in Fig. ^ 
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Figure 5: The scaled and normalised potential as a function of r/ro as obtained on 12^-24. 
The solid curve is the prediction from eq. (|^) . 


could be detected. Thus, there is at present no conclusive evidence for string breaking at 
length scales up to r ~ 1.5 fm. There are a number of arguments why this is so. Firstly, 
we wish to stress that the data points which probe the largest separations in Fig.|^ have 
typically been obtained using smaller values of for which the sea quarks may still 
be too heavy in order to produce a signihcantly different qualitative behaviour of V{r). 
It has also been argued @| that the Wilson loop used to extract V{r) does not project 
well onto states of broken strings. This has, in fact, been conhrmed in simulations using 
bosonic matter helds [^,^. In QCD, further investigations of this issue are required, in 
particular for smaller sea quark masses. Without a clear demonstration of string breaking, 
we can only give a rough estimate for the breaking scale Vb from the intersection of the 

From Fig. ^ we read off Vb — 2.6 tq. 


data of the potential with the value 2 




At small distances, where the potential is dominated by the Coulombic part, we hnd that 
the expression in eq. (P^D still describes the data surprisingly well, although the points 
obtained for the three lightest quark masses have a tendency to lie somewhat below the 
curve at the smallest separations. This implies that the data at small distances and 
quark masses seem to favour a larger value for e compared to e = 7r/12 in the pure gauge 
theory, as has been observed also in P|. This qualitative observation is consistent with 
the expected influence of dynamical quarks on the short-distance regime of the potential 
through the n/-dependence in the running of the strong coupling constant. 
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However, if one wants to quantify the change in e, one should take into account the lattice 
Greens function Giir) for one-gluon exchange in order to account for lattice artefacts at 
small distances. If tq/t in the last term of eq. (p^) is replaced by roGiir), it turns out 
that the deviation between the data points and the curve eq. for the two smallest 
values of r/ro is smaller but still signihcant. We estimate that the enhancement in e 
amounts to about 5-10% at the smallest value of the sea quark mass. 

This is only a crude analysis of sea quark effects in the short-distance part of the potential. 
In principle these effects on the running coupling could be probed by computing the 
coupling constant from the force according to 

«(^) = ^F{r)r^, (31) 

and comparing its scale dependence to the two-loop perturbative /3-function for nj = 2. 
In view of the many caveats concerning our present data, such as the fairly large length 
scales, the relatively heavy sea quarks and the lack of a continuum extrapolation, we have 
not seriously attempted such an analysis at this stage. 

To summarise, as far as the issue of string breaking is concerned, we hnd no hard evidence 
for the effects of dynamical quarks for distances up to 1.5 fm in the static quark potential. 
However, there are indications of a qualitatively different behaviour in the Coulombic 
range at small distances, which is hard to quantify and corroborate with the present data. 


4 Hadron spectroscopy 

In this section we describe the computation of the light hadron spectrum. The simulation 
parameters used have been discussed in subsect. p73[ 


4.1 Analysis and fitting procedure 


The amplitudes and masses of hadrons are obtained in a standard way by correlated least- 
fits of the correlation functions. The fitting function used was different for mesonic 
and baryonic channels. In the mesonic case, we have taken into account the backward 
propagating particle on a periodic lattice by fitting to the function 

Cuit) = -1- (32) 

where Ai,mi, i = 0,1 are the amplitudes and masses of the ground and first excited 
states. In the baryonic case, we have used 

Gsit) = (33) 


We have computed hadronic two-point correlation functions with different combinations 
of “fuzzing” both at source and at sink; we denote as FF the correlator fuzzed at 
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Figure 6: Effective mass plots on a 12^ ■ 24 volume for = 0.1398. In (a) 

the pseudoscalar and in (b) the nucleon is shown for different fuzzing combinations, i.e. 
FF (o), FL (O), LF (□) and LL (x). 


source and sink, FL fuzzed only at source, etc. We have found that the FF correlator 
allows the fastest isolation of the fundamental state, even in the case of the lightest k’s in 
which the effect of fuzzing is most important. To illustrate this point, we show in Fig. § 
the effective mass plots of the pseudoscalar and the nucleon for = 0.1398 on 

a 12^ ■ 24 volume. We conclude that the effect of fuzzing is quite signihcant compared 
with the LL case, as also found in the quenched approximation 10 . We have htted 
simultaneously the LL and FF correlators to a double exponential functional form, using 
the difference between the correlators to control the hrst excited state. 

The quantity y^/d.o.f. used to monitor the quality of a correlated ht is known to suffer 
from a systematic bias, which depends on the d.o.f. and the statistics H- We have imple¬ 
mented the technique of eigenvalue smoothing in the computation of the correlation 
matrix to take this bias into account. We have performed a “sliding window” analysis, 
hxing the maximum value of the ht interval, i.e. tmax = H, and varying to monitor the 
optimal interval. The stability criterion we have used is that the hadron mass should not 
change appreciably as tmin is changed by one unit. The value of tmin has been determined 
for each different combination of and k™*. 


4.2 The dynamical spectrum 

In a numerical simulation with dynamical fermions the parameters k®®'^ and are dis¬ 
tinct and each set of conhgurations generated for different is independent. We have 
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Figure 7: Effective mass plots for the pseudoscalar (a), vector (b), nucleon (c) and A (d), 
on a 12^ ■ 24 volume for all values of The effective mass has been calculated from the 
FF correlator with k®®®' = Solid lines denote the htted value of the mass, obtained 
from a simultaneous £t to the FF and LL correlators. The dashed lines denote the error 
band, and the length of the lines indicates the £t interval. 


performed the analysis of the hadron spectrum for each hxed value of At an in¬ 

termediate level, these simulations can be thought of as “pseudoquenched”, which come 
closer to the description of the real world as the sea quark mass approaches its physical 
value. There is another reason why we found it useful to simulate different values of 
for each cf. Table We interpret the heavier k’s as describing the valence quarks, 
in particular the strange, in the sea of light quarks, i.e. the up and the down. 

In Fig. 1^ we show the effective mass plots on a 12^ ■ 24 volume for the pseudoscalar, 
vector, nucleon and A, obtained for k®®®' = k™'. The pseudoscalar shows a clear plateau 
at all values of the hopping parameter, whereas for the vector the plateau becomes more 
unstable at the lightest quark masses. In the baryonic channels, on the other hand, the 
plateaux are, as expected, more fluctuating and we require a longer lattice in the time 
direction. 
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L = 

12 

L = 16 

t^sea 

hj 

mps/my 

mps/my 

0.1370 

0.855 

+ 2 
- 3 


0.1380 

0.825 

+ 4 
- 5 


0.1390 

0.785 

+ 4 
- 7 

0.785 + ® 

0.1395 

0.710 

+ 10 
-10 

0.719 

0.1398 

0.674 

+ 9 
-23 

0.670 +53 


Table 5: The ratio mps/my for = k™* on the two largest lattice sizes. 


In the tables given in Appendix ^ we summarise the results for the hadron masses ob¬ 
tained on the different volumes. We give the masses both in lattice units and in units of 
tq, the latter being more signihcant in the comparison between different volumes, as it 
compensates for the sea quark dependence of the lattice spacing. 

In Table D we list the values of the mass ratio of pseudoscalar and vector mesons, mps/my, 
obtained for which is a measure of how heavy our dynamical quark masses 

are relative to the real up and down quarks. Given that our lightest sea quark produces 
mps/my ~ 0.67, we conclude that the sea quarks used in our simulation are still relatively 
heavy. 


A useful quantity is the critical value of the hopping parameter, Kent- Here, we have used 
our data for pseudoscalar mesons, computed for and determined Kcrit from a 


£t 


(amps)^ = aB (airiq), 


(34) 


where 


mq = mq{l + b,e.amg), amg = ^(^- 

L \ rC ^crit / 


( 35 ) 

' ^ ^crit ^ 

Since 6m has not been determined non-perturbatively, we have used its perturbative ex¬ 
pression at one loop 

6m = -^- 0.0962^/2. (36) 


We considered this choice sufficient for our purposes, given that the value for Cs„ used in 
this study does not completely remove the leading cutoff effects. For every lattice size we 
fitted the three most chiral points to eq. (|3^), and the results for Kent are 


L/a = 

8 : 

Kent = 0.14047 

+ 6 
-7 

(37) 

L/a = 

12 : 

Kcrit = 0.14040 

+ 2 
-6 

(38) 

L/a = 

16 : 

«(crit = 0.14043 

+ 2 
- 3 ' 

(39) 
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As an aside we remark that the results for Kent obtained using 6m = 0 in eq. (|35| ) are 
entirely compatible with these values within errors. 

Traditionally, the way to make contact with the physical values of the light hadron spec¬ 
trum in quenched simulations has been to extrapolate the masses, obtained at several 
values of k, to the chiral limit. For example, the lattice spacing has usually been deter¬ 
mined by extrapolation of the vector mass and the physical value of nip. This approach is 
safe as long as the particles have zero decay width, as in the quenched approximation, but 
is no longer feasible in the dynamical case in which the p is not stable. Hence, it desirable 
to avoid extrapolations to the chiral limit, whenever possible. With this viewpoint, it has 
been proposed in |]4^, ^] to extract physical values from the region of the strange quark. 


4.3 Finite-size effects in the spectrum 


As already mentioned in the introduction, one of the aims of this study is to acquire 
experience of the systematics in dynamical simulations with improved fermions, even 
if the 0{a) effects are not entirely removed. One important feature to address is the 
presence of hnite-size effects in the spectrum, which determines the volume at which 
we can reliably carry out the calculation. Finite-size effects in dynamical simulations 
of hadronic spectrum with Wilson-like fermions have not been studied in great detail. 
The only results are those of the SESAM and TyL Collaborations with an unimproved 
fermionic action, exploring volumes 16^ ■ 32 and 24^ • 40 [|^ . It is thus important to study 
and quantify these effects using an 0(a)-improved action. 


In Figs. H and ^ we show the volume dependence of the meson and baryon masses, in 
units of ro, as a function of for It conhrms the behaviour already 

found in the study of the static quark potential. That is, we hnd pronounced hnite-size 
effects between 8^ ■ 24 and 12^ ■ 24, which can grow up to 15 — 17% in the mesonic sector 
and up to 25 — 28% in the baryonic sector, as we move towards the most chiral point. 
On the other hand, between 12^ ■ 24 and 16^ ■ 24 we hnd no signiheant discrepancy within 
statistical accuracy at all common values of the quark mass. From the similarity of the 
hnite-size behaviour of meson masses and the static quark potential discussed in section |^, 
we conclude that the bound eq. is also valid for simple hadronic quantities. 


4.4 Sea quark effects in the spectrum 


The parameter J has been introduced in [^] as a way to quantify the discrepancy between 
the quenched spectrum and experiment. It is dehned as 


dm^J 

J Hlx* — 2 : 


dm 


(40) 


PS 
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Figure 8: Vector and pseudoscalar masses, in units of rg, as a function of 
for different volumes. 




and has the attractive feature that it is dimensionless, does not involve extrapolations 
to the chiral limit and is independent of the quark mass values chosen to evaluate it, 
provided that my depends linearly on rupg. The only physical prediction sacrihced to 
calculate J is the ratio mx*/TnK = 1.81. In the (mpg,mv) plane this corresponds to a 
parabola, the intercept of which with the linear interpolation of the simulated data yields 
the point of strange mesons. An alternative way to determine J is to use 

mesons with full strange valence quark content, by assuming that the rjs is purely ss and 
that = 1.57. 


We emphasise that a realistic evaluation of J in the dynamical case is not straightforward. 
As pointed out in it would not be appropriate to compare meson masses obtained 
using different dynamical quark masses as the lattice spacing a depends on the sea quark 
mass. Our approach has been to £x the sea quark mass and for each sea quark consider 
different valence quark masses. Ideally, since we are looking to interpret the valence 
quarks as having strange flavour in the sea of light quarks, we would need to consider 
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Figure 9: A and nucleon masses, in units of ro, as a function of for different 

volumes. 


values of With our present dataset, this is only possible at the most chiral 

sea quark masses. However, since even our lightest sea quark mass is in the region of that 
of the strange quark, we do not expect the values of J to be signihcantly closer to the 
experimental value of J = 0.48(2), compared to the quenched approximation. 


In Fig. |T^ we show the values of J and compare them with the quenched result at /? = 5.7 
and V = 12^ ■ 24 |@. The values of J, obtained by either hxing mK^/mK or m^/rrir,^ 
show no appreciable trend towards the experimental point as the sea quark mass decreases. 
Moreover, the errors are amplihed by the htting process to extract the slope dmY/dmpg, 
especially at the lightest 


Another way to look for dynamical quark effects in the light hadron spectrum, which does 
not involve any htting, is to concentrate directly on the plot of the vector mass versus the 
pseudoscalar mass squared. This is shown in Fig. 0. As the sea quark mass is decreased 
relative to the valence quark mass, one observes a signihcant, albeit small trend of the 
meson masses towards the point i.e. the pair of mesons, whose valence quark 
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Figure 10: Values of the parameter J, determined from fixing the K*/K and plotted 
against for V = 12^ ■ 24. The quenched result, taken from [Q, is represented by 

the dashed lines. 


content resembles most closely that used in our simulation. Concentrating on the lighter 
sea-quark mass, e.g. = 0.1395, we can assert a signihcant shift compared to the 
quenched result, at /3 = 5.7 and V = 12^ ■ 24 p6 . 


At this stage it is hard to quantify the observed shift, and to disentangle the genuine sea 
quark effect from residual lattice artefacts, which could be fairly large in these simulations. 
A suitable approach would be to monitor dynamical quark effects at hxed lattice spacing. 
Starting from the quenched approximation, and going to ever lighter sea quark masses, 
one would have to perform a sequence of simulations, which are matched such that they all 
reproduce the same value of a suitable lattice scale, e.g. vq |^. Results will be published 
in a future publication, using the fully 0(a) improved action for dynamical and quenched 
simulations 


Another quantity, which can be used to highlight the effects of dynamical quarks, is 
the vector-pseudoscalar mass splitting. It is well known that lattice simulations fail to 
reproduce the experimental fact that this hyperhne splitting is constant over a large range 
of quark masses, rriy — rupg ~ 0.55 GeV^. The discrepancy between the experimental 
and (much lower) lattice determinations of rriy — rupg are partly due to lattice artefacts 
|48| . However, it is widely expected that the remaining difference, which becomes more 


pronounced as the valence quark mass is increased, is due to dynamical quark effects. 
In Fig. ^ we plot the meson hyperhne splitting (myro)^ — (mpsro)^ versus (mpsro)^ for 
all values of The numerical data are compared to the 4> — tJs splitting. Despite the 
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Figure 11: Vector mass plotted against the pseudoscalar mass for V = 12^ ■ 24 for 
eral sets of different sea quark masses. Asterisks denote the experimental points, 
quenched results (labelled = 0) are taken from PB|. 


sev- 

The 


relatively poor statistical accuracy it is obvious that the numerically determined hyperhne 
splitting shows a trend towards the experimental point as the sea quark mass is decreased. 
As regards the quantihcation of the relative deviation from experiment between the results 
obtained in the quenched approximation and for hnite sea quark mass, the same caveats 
apply as for the interpretation of Fig. |^, namely, that the comparison should be performed 
for hxed lattice spacing. 


5 Conclusions 

In this paper, we have presented the hrst results on the light hadron spectrum and 
the static quark potential obtained from dynamical simulations using an 0{a) improved 
fermion action at (3 = 5.2. Sea quark masses were chosen such that mps/my was in the 
range 0.86 down to 0.67. The value of = 1.76 was not appropriate to remove all 
leading discretisation effects. We wish to point out, however, that at such a low /3 value 
residual lattice artefacts could be relatively large, even after full 0(a) improvement. This 
question clearly requires further investigation. 

We have addressed the important issue of hnite-size effects, which are expected to be 
larger in dynamical simulations, compared to the quenched approximation. By simulating 
three different lattice sizes for a range of sea quark masses, we found that hnite-size 
effects are practically absent for box sizes L ^ 1.6 fm, and sea quark masses which give 
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Figure 12: Hyperfine splitting (myro)^—(rupsro)^ plotted against (mpsro)^ for V = 12^-24. 
The quenched results = 0) are taken from |^ . 

nT'Ps/fny ~ 0.67. This is observed in the data for the static quark potential and hadron 
masses. However, for lighter dynamical quarks one would expect that yet larger volumes 
are required. 

Instead of presenting quantitative results for the light hadron spectrum, we have concen¬ 
trated on highlighting the effects of dynamical quarks. Although the evidence for string 
breaking in the data for the static quark potential and for an improved behaviour of the 
parameter J remains inconclusive, we have detected significant effects due to dynamical 
quarks. Namely, the Coulombic part of the static quark potential is enhanced for finite 
sea quark mass, and the vector-pseudoscalar hyperfine splitting moves closer to the exper¬ 
imental value as the sea quark mass is decreased. Furthermore, pairs of (myro, (rupsro)^) 
show a trend towards experiment when dynamical quarks are “switched on”. 

The results presented here serve as a guideline for ongoing investigations, which are per¬ 
formed using the fully 0(a) improved action for two flavours of dynamical quarks, whose 
masses are closer to the chiral limit than those used in this study. 
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A Tables 


Here we list the meson and baryon masses for all combinations of and k™*, and on 
all lattice sizes used in the simulations. 

In the following tables we also list the htting interval and the value of x^/d.o.f., obtained 
from a correlated £t to eqs. (|3^) or (^), respectively. Some of the hts produce values of 
X^/d.o.f., which would normally be regarded as unacceptably large, even though we are 
conhdent that the £t intervals were properly tuned. In such cases we have compared the 
results for the masses to those obtained using an uncorrelated £t over the same interval. 
We always found that the uncorrelated £t gave results which, within errors, were perfectly 
compatible with those from the correlated £t, whilst producing values for y^/d.o.f., which 
were signihcantly below one. Thus we are conhdent that the results presented in the tables 
below are reliable, that a double exponential is an appropriate model function in all cases, 
but that correlations in the data are clearly present. 


L3 

■T 

i^sea 

hj 

^val 

rt 

mpsa 

mpsro 

Fit 

xV 

dof 

8^ 

24 

0.1370 

0.1370 

1.119 

+ 8 
- 4 

2.502 

+47 

-52 

[ 6.11] 

12.70 

/6 




0.1380 

1.059 

+ 8 
- 4 

2.368 

+45 

-50 

[ 6,11) 

12.77 

/6 




0.1390 

0.997 

+ 9 
- 5 

2.229 

+44 

-47 

1 6.111 

13.20 

/6 




0.1395 

0.965 

+ 9 
- 4 

2.158 

+43 

-45 

[ 6.11] 

13.57 

/6 

8^ 

24 

0.1380 

0.1370 

1.011 

+ 10 
- 2 

2.502 

+79 

-68 

[ 5.11] 

14.94 

/8 




0.1380 

0.945 

+ 11 

- 2 

2.339 

+75 

-64 

[ 5.111 

13.28 

/8 




0.1390 

0.874 

+ 12 
- 6 

2.164 

+72 

-61 

[ 5.11] 

7.26 

/8 




0.1395 

0.840 

+ 14 

- 6 

2.078 

+71 

-58 

1 5.111 

7.50 

/8 

8^ 

24 

0.1390 

0.1370 

0.870 

+ 9 
- 5 

2.514 

+58 

-53 

[ 6.11] 

8.14 

/6 




0.1380 

0.799 

+ 10 
- 5 

2.309 

+57 

-49 

1 6.111 

8.18 

/6 




0.1390 

0.726 

+ 11 
- 6 

2.099 

+53 

-45 

[ 5.111 

10.24 

/8 




0.1395 

0.686 

+ 12 
- 5 

1.982 

+53 

-43 

[ 5.11) 

10.86 

/8 

8^ 

24 

0.1395 

0.1370 

0.786 

+ 9 
-14 

2.924 

+78 

-78 

[ 6.11) 

9.70 

/6 




0.1380 

0.710 

+ 9 
-12 

2.640 

+72 

-68 

[ 6.11] 

10.75 

/6 




0.1390 

0.642 

+ 15 
-10 

2.387 

+80 

-60 

[ 5.111 

15.07 

/8 




0.1395 

0.603 

+ 17 
- 7 

2.242 

+82 

-51 

[ 4.111 

17.30 

/lO 


Table 6: Pseudoscalar masses on 8^ ■ 24, in lattice units and in units of tq. 
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CO 

^sea 

hj 

^val 

hj 

mpsa 

mpsro 

Fit 

xV dof 

12^ ■ 24 

0.1370 

0.1370 

1.108 

+ 3 

- 2 

2.541 

+23 

-26 

[6,11] 

9.93 / 6 



0.1380 

1.047 

+ 3 

- 2 

2.402 

+22 

-25 

15,11] 

11.18 / 8 



0.1390 

0.983 

+ 4 

- 2 

2.255 

+21 

-23 

|5,11] 

11.65 / 8 



0.1395 

0.950 

+ 4 
- 3 

2.179 

+21 

-23 

|5,11] 

11.83 / 8 

123.24 

0.1380 

0.1370 

1.000 

+ 3 
- 2 

2.569 

+35 

-30 

[6,11] 

10.41 / 6 



0.1380 

0.936 

+ 3 
- 2 

2.404 

+33 

-29 

[ 6,11] 

10.79 / 6 



0.1390 

0.869 

+ 4 

- 2 

2.232 

+31 

-27 

[6,11] 

11.37 / 6 



0.1395 

0.834 

+ 4 

- 2 

2.142 

+30 

-26 

|6,11] 

11.56 / 6 

123.24 

0.1390 

0.1370 

0.858 

+ 4 

- 3 

2.613 

+51 

-45 

[ 7,11] 

5.93 / 4 



0.1380 

0.785 

+ 4 

- 3 

2.391 

+48 

-42 

|7,11] 

4.16 / 4 



0.1390 

0.707 

+ 5 
- 3 

2.155 

+44 

-38 

|7,11] 

2.69 / 4 



0.1395 

0.669 

+ 5 
- 3 

2.037 

+42 

-36 

[ 6,11] 

16.86 / 6 

123.24 

0.1395 

0.1370 

0.775 

+ 6 
- 4 

2.662 

+42 

-39 

[ 7,11] 

21.46 / 4 



0.1380 

0.699 

+ 5 

- 4 

2.402 

+38 

-36 

|6,11] 

26.14 / 6 



0.1390 

0.615 

+ 6 
- 4 

2.112 

+36 

-32 

[ 6,11] 

19.70 / 6 



0.1395 

0.558 

+ 8 
- 7 

1.916 

+38 

-35 

|5,11] 

8.96 / 8 

123.24 

0.1398 

0.1380 

0.639 

+ 5 
- 4 

2.332 

+27 

-21 

[ 7,11] 

3.63 / 4 



0.1390 

0.551 

+ 4 

- 9 

2.011 

+23 

-37 

|7,11] 

2.09 / 4 



0.1395 

0.492 

+ 3 
- 9 

1.797 

+ 18 
-36 

|6,11] 

3.41 / 6 



0.1398 

0.476 

+ 3 
-16 

1.738 

+ 17 
-61 

[6,11] 

8.14 / 6 


Table 7: Pseudoscalar masses on 12^ ■ 24, in lattice units and in units of tq. 












CO 

^sea 

hj 

i^val 

hj 

mpsa 

mpsro 

Fit 

xV dof 

16^ ■ 24 

0.1390 

0.1390 

0.701 

— 5 

2.120 

+30 

-22 

[6,11] 

10.62 / 6 



0.1395 

0.660 +4 

1.998 

+28 

-21 

[6,11] 

9.31 / 6 



0.1398 

0.635 

— 5 

1.922 

+27 

-22 

[6,11] 

8.34 / 6 

16^ ■ 24 

0.1395 

0.1390 

0.610 +2 

2.101 

+28 

-36 

[6,11] 

7.04 / 6 



0.1395 

0.564 

1.942 

+25 

-36 

[6,11] 

3.05 / 6 



0.1398 

0.537 +5 

1.848 

+26 

-32 

[7.11] 

2.99 / 4 

16^ ■ 24 

0.1398 

0.1390 

0.551 

— 5 

2.011 

+26 

-26 

[ 7,11] 

10.66 / 4 



0.1395 

0.502 

1.834 

+24 

-28 

17,11] 

8.48 / 4 



0.1398 

0.468 

— 5 

1.707 

+23 

-24 

[7,11] 

6.02 / 4 

ble 8: Pseudoscalar masses 

on 16^■24 

in lattice units and 

in units of 

L^-T 

i^sea 

hj 

^val 

hj 

mva 

myro 

Fit 

xV dof 

8^ ■ 24 

0.1370 

0.1370 

1.305 

— D 

2.917 

+ 57 
- 61 

] 7,11] 

5.82 / 4 



0.1380 

1.257 

— D 

2.810 

+ 56 
- 59 

] 7,11] 

5.09 / 4 



0.1390 

1.209 

2.703 

+ 56 
- 58 

] 7,11] 

5.27 / 4 



0.1395 

1.185 

2.650 

+ 56 
- 57 

] 7,11] 

5.78 / 4 

8^ ■ 24 

0.1380 

0.1370 

1.190 

2.946 

+ 95 
- 82 

] 6,11] 

7.30 / 6 



0.1380 

1.140 

2.822 

+ 92 
- 79 

] 6,11] 

7.68 / 6 



0.1390 

1.090 

2.697 

+ 90 
- 76 

[6,11] 

9.02 / 6 



0.1395 

1.064 

2.634 

+ 89 
- 76 

] 6,11] 

10.05 / 6 

8^ ■ 24 

0.1390 

0.1370 

1.045 

3.022 

+ 72 
- 64 

] 7,11] 

5.21 / 4 



0.1380 

0.998 

2.885 

+ 76 
- 64 

1 6,11] 

6.49 / 6 



0.1390 

0.944 

2.730 

+ 79 
- 61 

] 6,11] 

4.66 / 6 



0.1395 

0.918 

2.653 

+ 80 
- 61 

]6,11] 

4.08 / 6 

8^ ■ 24 

0.1395 

0.1370 

0.961 +^g 

3.573 

+141 
- 77 

] 7,11] 

12.79 / 4 



0.1380 

0.917 

3.411 

+122 
- 70 

] 6,11] 

11.80 / 6 



0.1390 

0.874 

3.251 

+149 

- 65 

] 6,11] 

11.01 / 6 



0.1395 

0.850 +^g 

3.159 

+124 
- 69 

] 6,11] 

11.17 / 6 


Table 9: Vector masses on 8^ ■ 24, in lattice units and in units of tq. 



















CO 

^sea 

hj 

^val 

hj 

my a 

mv+o 

Fit 

xV dof 

12^ ■ 24 

0.1370 

0.1370 

1.296 

+ 6 

- 4 

2.973 

+29 

-31 

[ 7.11] 

3.85 / 4 



0.1380 

1.249 

+ 6 
- 5 

2.866 

+28 

-31 

1 7,11] 

4.17 / 4 



0.1390 

1.202 

+ 6 
- 5 

2.757 

+28 

-30 

|7,11] 

4.77 / 4 



0.1395 

1.178 

+ 7 
- 5 

2.703 

+28 

-29 

1 7,11] 

5.19 / 4 

123.24 

0.1380 

0.1370 

1.185 

+ 6 
- 4 

3.043 

+43 

-37 

[ 7,11] 

24.88 / 4 



0.1380 

1.135 

+ 7 
- 4 

2.915 

+43 

-36 

[ 6,11] 

33.99 / 6 



0.1390 

1.084 

+ 8 
- 4 

2.783 

+42 

-34 

[6,11] 

28.64 / 6 



0.1395 

1.058 

+ 9 
- 4 

2.717 

+42 

-34 

16.11] 

25.81 / 6 

123.24 

0.1390 

0.1370 

1.016 

+ 4 

- 5 

3.094 

+60 

-55 

[6,11] 

6.99 / 6 



0.1380 

0.962 

+ 7 
- 3 

2.931 

+59 

-51 

|7,11] 

14.24 / 4 



0.1390 

0.901 

+ 10 
- 3 

2.746 

+61 

-48 

|7,11] 

8.29 / 4 



0.1395 

0.891 

+ 9 
- 1 

2.714 

+59 

-47 

[ 6,11] 

36.73 / 6 

123.24 

0.1395 

0.1370 

0.930 

+ 9 
- 7 

3.195 

+54 

-49 

[ 7,11] 

29.92 / 4 



0.1380 

0.867 

+ 7 
- 6 

2.980 

+48 

-46 

|6,11] 

12.76 / 6 



0.1390 

0.814 

+ 8 
- 6 

2.796 

+48 

-43 

|5,11] 

13.67 / 8 



0.1395 

0.786 

+ 9 
- 7 

2.700 

+49 

-44 

|5,11] 

13.10 / 8 

123.24 

0.1398 

0.1380 

0.791 

+ 16 
- 1 

2.890 

+64 

-20 

[ 7,11] 

25.71 / 4 



0.1390 

0.735 

+ 7 

-14 

2.685 

+33 

-54 

|6,11] 

14.22 / 6 



0.1395 

0.725 

+ 12 

-16 

2.647 

+49 

-62 

|5,11] 

29.11 / 8 



0.1398 

0.706 

+ 12 

-16 

2.578 

+47 

-61 

|5,11] 

25.26 / 8 


Table 10: Vector masses on 12^ ■ 24, in lattice units and in units of tq. 













T 

hj 

^val 

hj 

mya 

mv+o 

Fit 

xV dof 

16^ 

24 

0.1390 

0.1390 

0.893 

+ 9 
- 6 

2.702 

+39 

-28 

[ 7.11] 

14.64 / 4 




0.1395 

0.863 

+ 10 
- 6 

2.610 

+40 

-28 

1 7,11| 

14.79 / 4 




0.1398 

0.852 

+ 11 
- 6 

2.580 

+43 

-27 

[4,11] 

22.82 / 10 

16^ 

24 

0.1395 

0.1390 

0.813 

+ 8 
- 6 

2.801 

+43 

-51 

[ 7,11] 

13.43 / 4 




0.1395 

0.785 

+ 9 
- 7 

2.702 

+45 

-51 

[ 7,11| 

12.70 / 4 




0.1398 

0.773 

+ 10 
- 8 

2.664 

+46 

-53 

[ 5.11] 

14.00 / 8 

16^ 

24 

0.1398 

0.1390 

0.749 

+ 9 
- 8 

2.733 

+41 

-37 

1 7,11| 

2.09 / 4 




0.1395 

0.716 

+ 10 
- 8 

2.614 

+44 

-37 

1 7,11] 

1.30 / 4 




0.1398 

0.698 

+ 11 
-10 

2.547 

+47 

-41 

[ 7,11] 

1.60 / 4 


Table 11: Vector masses on 16^ ■ 24, in lattice units and in units of tq. 


L'i . J' I^sea I^val friNCl mArro Fit X^/ dof 


8^ ■ 24 

0.1370 

0.1370 

2.094 

+42 

-14 

4.681 

+124 

-102 

1 6.11] 

14.59 / 6 



0.1380 

2.020 

+43 

-15 

4.517 

+124 
- 99 

[ 6,11] 

13.29 / 6 



0.1390 

1.945 

+44 

-18 

4.349 

+124 
- 98 

[ 6,11] 

12.07 / 6 



0.1395 

1.907 

+49 

-21 

4.264 

+132 
- 99 

[ 6,11] 

11.54 / 6 

8^ ■ 24 

0.1380 

0.1370 

1.926 

+40 

-14 

4.768 

+174 

-134 

[ 6,11] 

11.13 / 6 



0.1380 

1.830 

+47 

-10 

4.528 

+178 

-125 

[ 6,11] 

9.43 / 6 



0.1390 

1.730 

+55 
- 8 

4.282 

+188 

-117 

1 6,11] 

7.65 / 6 



0.1395 

1.677 

+63 
- 7 

4.151 

+198 

-114 

[ 6,11] 

6.74 / 6 

8^ ■ 24 

0.1390 

0.1370 

1.687 

+29 

-12 

4.877 

+132 

-103 

[ 6,11] 

11.89 / 6 



0.1380 

1.596 

+38 

-10 

4.613 

+146 
- 97 

1 6,11] 

9.77 / 6 



0.1390 

1.511 

+44 

-14 

4.368 

+155 
- 96 

[ 6,11] 

8.65 / 6 



0.1395 

1.473 

+47 

-20 

4.259 

+162 

-103 

[ 6,11] 

8.89 / 6 

8^ ■ 24 

0.1395 

0.1370 

1.606 

+25 

-22 

5.971 

+171 

-141 

[ 6,11] 

8.73 / 6 



0.1380 

1.535 

+24 

-26 

5.706 

+164 

-148 

[ 6,11] 

7.84 / 6 



0.1390 

1.463 

+34 

-48 

5.440 

+181 

-209 

[ 7,11] 

4.81 / 4 



0.1395 

1.397 

+42 

-26 

5.195 

+200 

-140 

[4,11] 

9.34 / 100 


Table 12: Nucleon masses on 8^ ■ 24, in lattice units and in units of tq. 



















L3 

T 

^sea 

hj 

^val 

hj 

7TIn(1 

niNVo 

Fit 

xV dof 

12^ 

24 

0.1370 

0.1370 

2.053 

+32 

-22 

4.710 

+ 85 
- 69 

1 7.11] 

8.29 / 4 




0.1380 

1.975 

+36 

-19 

4.531 

+ 91 
- 62 

[ 7,111 

7.97 / 4 




0.1390 

1.900 

+33 

-10 

4.359 

+ 85 
- 50 

[6,11] 

10.42 / 6 




0.1395 

1.860 

+33 

-11 

4.267 

+ 84 
- 49 

[6,111 

9.47 / 6 

123 

24 

0.1380 

0.1370 

1.850 

+24 
- 9 

4.752 

+ 88 
- 60 

1 7,111 

4.95 / 4 




0.1380 

1.763 

+23 

-10 

4.528 

+ 85 
- 58 

[ 7,111 

4.85 / 4 




0.1390 

1.673 

+24 
- 9 

4.296 

+ 83 
- 56 

1 7,111 

5.56 / 4 




0.1395 

1.626 

+24 

-11 

4.176 

+ 83 
- 56 

[ 7,111 

6.19 / 4 

123 

24 

0.1390 

0.1370 

1.604 

+20 
- 9 

4.886 

+111 

- 88 

1 7,111 

15.34 / 4 




0.1380 

1.505 

+21 

-10 

4.585 

+107 
- 83 

[ 7,111 

13.01 / 4 




0.1390 

1.399 

+22 

-11 

4.263 

+106 
- 80 

1 7,111 

10.61 / 4 




0.1395 

1.343 

+24 

-11 

4.089 

+107 
- 77 

1 7,111 

9.33 / 4 

123 

24 

0.1395 

0.1370 

1.434 

+ 19 
-10 

4.927 

+ 95 
- 76 

1 7,111 

13.69 / 4 




0.1380 

1.335 

+20 

-11 

4.585 

+ 94 
- 73 

[ 7,111 

10.29 / 4 




0.1390 

1.233 

+20 

-12 

4.234 

+ 92 
- 71 

[ 7,111 

7.69 / 4 




0.1395 

1.178 

+22 

-10 

4.045 

+ 95 
- 65 

[ 6,11] 

5.05 / 6 

123 

24 

0.1398 

0.1380 

1.236 

+26 

-15 

4.514 

+102 
- 63 

|6,11| 

15.02 / 6 




0.1390 

1.124 

+23 

-23 

4.106 

+ 90 
- 89 

[ 7,111 

14.30 / 4 




0.1395 

1.069 

+25 

-20 

3.902 

+ 97 
- 76 

[6,111 

14.69 / 6 




0.1398 

1.031 

+26 

-20 

3.765 

+101 

- 78 

[ 6,11] 

14.74 / 6 


Table 13: Nucleon masses on 12^ ■ 24, in lattice units and in units of tq. 












^3 . rj. ^sea 




UInTo 


K 


val 


Fit 


xV dof 


16^ 

■ 24 

0.1390 

0.1390 

1.345 

+33 
- 5 

4.071 

+110 
- 36 

[ 7.11] 

5.49 

/4 




0.1395 

1.295 

+32 
- 6 

3.918 

+106 
- 36 

1 7,11] 

6.24 

/4 




0.1398 

1.264 

+31 
- 7 

3.826 

+102 
- 37 

1 7,11] 

6.72 

/4 

16^ 

■ 24 

0.1395 

0.1390 

1.261 

+25 

-10 

4.344 

+ 99 
- 80 

[ 7,11] 

13.78 

/4 




0.1395 

1.205 

+24 

-11 

4.149 

+ 96 
- 79 

1 7,11] 

10.92 

/4 




0.1398 

1.164 

+24 

-19 

4.010 

+ 96 
- 92 

1 7,11] 

9.08 

/4 

16^ 

■24 

0.1398 

0.1390 

1.144 

+36 

-31 

4.178 

+136 

-117 

[ 7,11] 

23.51 

/4 




0.1395 

1.065 

+47 

-35 

3.888 

+173 

-133 

1 7,11] 

19.78 

/4 




0.1398 

1.030 

+43 

-34 

3.762 

+160 

-127 

1 7,11] 

17.48 

/4 


Table 14: 

Nucleon 

masses 

on 16^ 

24, 

in lattice units and in 

units of rx 

L^-T 

^sea 

hj 

t^val 

hj 

niAa 

niAro 

Fit 

xV dof 

8^ ■ 24 

0.1370 

0.1370 

2.146 

+39 

-13 

4.799 

+121 

-103 

[ 7.11] 

5.73 / 4 



0.1380 

2.070 

+40 

-14 

4.629 

+120 

-100 

[ 7,11] 

4.79 / 4 



0.1390 

1.989 

+44 

-13 

4.448 

+125 
- 96 

[ 7,11] 

4.64 / 4 



0.1395 

1.947 

+46 

-12 

4.353 

+128 
- 94 

[ 7,11] 

4.75 / 4 

8^ ■ 24 

0.1380 

0.1370 

2.131 

+51 
- 1 

5.274 

+202 

-143 

[ 6.11] 

16.72 / 6 



0.1380 

1.918 

+37 
- 4 

4.747 

+169 

-129 

[ 7,11] 

6.51 / 4 



0.1390 

1.854 

+34 
- 8 

4.588 

+160 

-126 

[ 7,11] 

3.99 / 4 



0.1395 

1.831 

+33 

-12 

4.531 

+158 

-126 

[ 7,11] 

3.93 / 4 

8^ ■ 24 

0.1390 

0.1370 

1.790 

+34 
- 9 

5.174 

+145 

-107 

[ 6,11] 

17.07 / 6 



0.1380 

1.717 

+32 
- 8 

4.964 

+139 

-102 

[ 6,11] 

14.01 / 6 



0.1390 

1.651 

+27 
- 9 

4.773 

+127 

-100 

[ 6,11] 

11.26 / 6 



0.1395 

1.621 

+27 

-13 

4.687 

+124 

-101 

[ 6,11] 

10.62 / 6 

8^ ■ 24 

0.1395 

0.1370 

1.718 

+42 
- 6 

6.389 

+219 

-126 

[ 6,11] 

10.07 / 6 



0.1380 

1.644 

+42 
- 6 

6.111 

+214 

-121 

[ 6,11] 

9.26 / 6 



0.1390 

1.578 

+38 

-13 

5.866 

+200 

-124 

[ 6,11] 

10.38 / 6 



0.1395 

1.538 

+43 

-19 

5.718 

+210 

-130 

[ 7,11] 

7.59 / 4 


Table 15: A masses on 8^ ■ 24, in lattice units and in units of tq. 



















L3. 

T 

^sea 

hj 

^val 

hj 

m/^a 

niAro 

Fit 

xV< 

iof 

12^ . 

■ 24 

0.1370 

0.1370 

2.139 

+29 

- 1 

4.907 

+ 79 
- 49 

[ 7.11] 

10.07 

/ 

4 




0.1380 

2.071 

+27 
- 3 

4.750 

+ 75 
- 48 

1 7,11] 

10.48 

/ 

4 




0.1390 

2.000 

+25 
- 4 

4.587 

+ 70 
- 47 

1 7,11] 

10.97 

/ 

4 




0.1395 

1.966 

+23 
- 5 

4.511 

+ 65 
- 47 

1 7,11] 

11.01 

/ 

4 

123 , 

■ 24 

0.1380 

0.1370 

1.933 

+30 
- 1 

4.963 

+101 
- 58 

[ 7,11] 

8.07 

/ 

4 




0.1380 

1.855 

+30 
- 2 

4.763 

+ 99 
- 56 

1 7,11] 

6.26 

/ 

4 




0.1390 

1.773 

+32 
- 2 

4.553 

+102 
- 53 

1 7,11] 

5.89 

/ 

4 




0.1395 

1.754 

+22 

-14 

4.504 

+ 82 
- 64 

1 7,11] 

5.19 

/ 

4 

123 , 

■ 24 

0.1390 

0.1370 

1.714 

+31 
- 5 

5.221 

+136 
- 90 

[ 7,11] 

18.87 

/ 

4 




0.1380 

1.630 

+30 
- 6 

4.964 

+131 

- 87 

1 7,11] 

16.59 

/ 

4 




0.1390 

1.545 

+27 
- 6 

4.707 

+123 
- 83 

1 7,11] 

13.88 

/ 

4 




0.1395 

1.502 

+29 
- 6 

4.576 

+124 
- 80 

1 7,11] 

12.39 

/ 

4 

123 , 

■ 24 

0.1395 

0.1370 

1.539 

+26 

-10 

5.285 

+115 
- 80 

[ 7,11] 

26.22 

/ 

4 




0.1380 

1.457 

+20 

-16 

5.005 

+ 97 
- 87 

1 7,11] 

21.69 

/ 

4 




0.1390 

1.372 

+ 19 
-19 

4.711 

+ 92 
- 91 

1 7,11] 

17.49 

/ 

4 




0.1395 

1.329 

+22 

-20 

4.565 

+ 99 
- 93 

1 7,11] 

15.43 

/ 

4 

123 , 

■ 24 

0.1398 

0.1380 

1.373 

+34 

-21 

5.016 

+129 
- 83 

[ 7,11] 

34.34 

/ 

4 




0.1390 

1.273 

+32 

-26 

4.648 

+121 

- 99 

1 7,11] 

27.18 

/ 

4 




0.1395 

1.222 

+32 

-27 

4.463 

+121 

-103 

1 7,11] 

23.00 

/ 

4 




0.1398 

1.191 

+29 

-28 

4.349 

+113 

-105 

1 7,11] 

20.55 

/ 

4 


Table 16; A masses on 12^ ■ 24, in lattice units and in units of tq. 












L^-T 

^sea 

^val 

hj 

niAa 

16^ ■ 24 

0.1390 

0.1390 

1.599 

+34 

-20 



0.1395 

1.555 

+34 

-19 



0.1398 

1.532 

+35 

-19 

16^ ■ 24 

0.1395 

0.1390 

1.395 

+28 
- 5 



0.1395 

1.347 

+30 
- 6 



0.1398 

1.318 

+30 
- 8 

16^ ■ 24 

0.1398 

0.1390 

1.296 

+29 

-11 



0.1395 

1.244 

+27 

-11 



0.1398 

1.208 

+29 

-12 


Table 17; A masses on 16^ ■ 24, in 


m/\rQ Fit / dof 
4.838 +^3 [ 6,11] 32.43 / 6 

4.704 [ 6,11] 28.16 / 6 

4.635 [ 6,11] 25.28 / 6 

4.805 [ 7,11] 11.43 / 4 

4.640 [ 7,11] 10.63 / 4 

4.538 [ 7,11] 10.20 / 4 

4.733 +^3 [ 7,11] 19.39 / 4 

4.541 [ 7,11] 17.22 / 4 

4.411 [ 7,11] 15.52 / 4 


units and in units of Tq. 













